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Abstract 

A relativistic nuclear energy density functional is developed, guided by two im- 
portant features that establish connections with chiral dynamics and the symmetry 
breaking pattern of low-energy QCD: a) strong scalar and vector fields related to in- 
medium changes of QCD vacuum condensates; b) the long- and intermediate-range 
interactions generated by one-and two-pion exchange, derived from in-medium chi- 
ral perturbation theory, with explicit inclusion of A(1232) excitations. Applications 
are presented for binding energies, radii of proton and neutron distributions and 
other observables over a wide range of spherical and deformed nuclei from 16 O to 
210 Po. Isotopic chains of Sn and Pb nuclei are studied as test cases for the isospin 
dependence of the underlying interactions. The results are at the same level of quan- 
titative comparison with data as the best phenomenological relativistic mean-field 
models. 
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1 Introduction. 



The successes of modern nuclear structure models in predicting many new 
phenomena in regions of exotic nuclei far from stability, and the recent appli- 
cations of chiral effective field theory to nucleon-nucieon scattering and the 
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few-body problem, have highlighted one of the fundamental problems in theo- 
retical nuclear physics: the relationship between low-energy, non-perturbative 
QCD and the rich structure of nuclear many-body systems. 

The most complete and accurate description of structure phenomena in heavy 
nuclei is currently provided by self-consistent non- relativist ic and relativistic 
mean-field approaches. By employing phenomenological effective interactions, 
adjusted to empirical properties of symmetric and asymmetric nuclear matter, 
and to bulk properties of stable nuclei, self-consistent mean-field methods have 
achieved a high level of precision in describing ground states and properties of 
excited states in stable nuclei, exotic nuclei far from /^-stability, and in nuclear 
systems at the nucleon drip-lines [1,2]. 

The self-consistent mean-field approach to nuclear structure represents an ap- 
proximate implementation of Kohn-Sham density functional theory (DFT) 
[3,4,5,6]. The DFT provides a description of the nuclear many-body problem 
in terms of an energy density functional, E[p\. Mean-field models approxi- 
mate the exact energy functional which includes all higher-order correlations. 
A major goal of nuclear structure theory is to build an energy density func- 
tional which is universal [3], in the sense that the same functional is used 
for all nuclei, with the same set of parameters. This framework should then 
provide a reliable microscopic description of infinite nuclear and neutron mat- 
ter, ground-state properties of bound nuclei, rotational spectra, low-energy 
vibrations and large-amplitude adiabatic properties [7]. 

In order to formulate a microscopic nuclear energy density functional, one 
must be able to go beyond the mean-field approximation and systematically 
calculate the exchange-correlation part, E xc [p\, of the energy functional, start- 
ing from the relevant active degrees of freedom at low energy. The exact E xc 
includes all many-body effects. Thus the usefulness of DFT crucially depends 
on our ability to construct accurate approximations to the exact exchange- 
correlation energy. The natural microscopic framework is chiral effective field 
theory. It is based on the separation of scales between long-range pion-nucleon 
dynamics, described explicitly, and short-distance interactions not resolved in 
detail at low energies. 

An extensive program, synthesizing effective field theory methods and density 
functional theory, has recently been initiated by Furnstahl and collaborators 
[8,9,10]. For a dilute, confined Fermi system with short-range interactions, an 
effective action formalism leads to a Kohn-Sham density functional by apply- 
ing an inversion method order-by-order in the relevant small expansion pa- 
rameter, the local Fermi momentum times the scattering length. The starting 
point is a generating functional with external sources coupled to composite 
Fermion operators, e.g. particle number densities and kinetic energy densi- 
ties. A functional Legendre transformation with respect to source fields leads 
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to an effective action functional from which the energy density functional is 
calculated. 

An alternative approach to the nuclear energy density functional, emphasizing 
links with low-energy QCD and its symmetry breaking pattern, has recently 
been introduced [11,12]. It is based on the following conjectures: 

(1) The nuclear ground state is characterized by strong scalar and vector 
mean fields which have their origin in the in-medium changes of the scalar 
quark condensate (the chiral condensate) and of the quark density. 

(2) Nuclear binding and saturation arise primarily from chiral (pionic) fluc- 
tuations (reminiscent of van der Waals forces) in combination with Pauli 
blocking effects and three-nucleon (3N) interactions, superimposed on 
the condensate background fields and calculated according to the rules 
of in-medium chiral perturbation theory (ChPT). 

The starting point is the description of nuclear matter based on the chiral ef- 
fective Lagrangian with pions and nucleons [13,14,15,16], recently improved by 
including explicit A(1232) degrees of freedom [17]. The relevant "small" scales 
are the Fermi momentum kf, the pion mass m n and the A — iV mass difference 
A = Ma — Mjv — 2.1m n , all of which are well separated from the characteris- 
tic scale of spontaneous chiral symmetry breaking, 4irf n ~ 1.16 GeV with the 
pion decay constant f n = 92.4 MeV. The calculations have been performed 
to three-loop order in the energy density. They incorporate the one-pion ex- 
change Fock term, iterated one-pion exchange and irreducible two-pion ex- 
change, including one or two intermediate A's. The resulting nuclear matter 
equation of state is given as an expansion in powers of the Fermi momentum 
kf. The expansion coefficients are functions of kf/m^ and A/m n , the dimen- 
sionless ratios of the relevant small scales. Divergent momentum space loop 
integrals are regularized by introducing subtraction constants in the spectral 
representations of these terms [17]. The (few) subtraction constants are the 
only parameters in this approach. They equivalently correspond to two- and 
three-nucleon contact interactions (and derivatives thereof), encoding short- 
distance dynamics not resolved in detail at the characteristic momentum scale 
kf <C 4:iT f n . The finite parts of the energy density, written in closed form as 
functions of kf/m n and A/m n , represent long and intermediate range (chiral) 
dynamics with input fixed entirely in the irN sector. The low-energy constants 
(contact terms) are adjusted to reproduce basic properties of symmetric and 
asymmetric nuclear matter. 

A first version (not yet including explicitly the A(1232)) of this microscopic 
approach has been tested in the analysis of bulk and single-nucleon properties 
of finite nuclei [11,12]. It was shown that chiral (two-pion exchange) fluctu- 
ations play a prominent role in nuclear binding and in the saturation mech- 
anism, while additional strong scalar and vector mean fields of about equal 



3 



magnitude and opposite sign, induced by changes of the QCD vacuum in the 
presence of baryonic matter, generate the large effective spin-orbit potential in 
finite nuclei. A first series of promising results for N m Z nuclei demonstrated 
that such an approach to nuclear dynamics, constrained by the chiral sym- 
metry breaking pattern and the condensate structure of low-energy QCD, has 
the capability of describing finite nuclei and their properties at a quantitative 
level comparable with phenomenological self-consistent mean-field models. 

Chiral effective field theories are systematically improved by introducing ex- 
plicit A(1232) degrees of freedom. Much better isospin properties of nuclear 
matter are found by including chiral 7riVA-couplings [17]. This has an amelio- 
rating influence on the isovector channel of the nuclear energy density func- 
tional for finite nuclei, and much improved results are expected for ground- 
state properties of N ^ Z nuclei. In the present work the effects of the inclu- 
sion of chiral 7riVA-dynamics on the nuclear energy density functional will be 
investigated. Specifically, the chiral nuclear matter energy density functional 
will be mapped onto the exchange-correlation energy density functional of a 
covariant point-coupling model for finite nuclei, including gradient corrections. 
This model will be employed in the description of ground-state properties of a 
broad range of spherical and deformed nuclei. The results will be analyzed in 
comparison with experimental data on binding energies, charge radii, neutron 
radii and deformation parameters for several isotopic chains. 

In Section 2 we construct the nuclear energy density functional based on the 
conjectures mentioned previously. Next, in order to deal with a broad range of 
finite (medium-heavy and heavy) nuclei, it is convenient to formulate an equiv- 
alent covariant point-coupling model with density-dependent contact interac- 
tions. The mapping of the nucleon self-energies in nuclear matter, calculated 
using chiral dynamics, onto those of the point coupling model for finite nu- 
clei and the fine-tuning of the remaining parameters is described in Section 3. 
In Section 4 the resulting self-consistent equations are solved for ground-state 
properties of a number of spherical and deformed nuclei. Section 5 summarizes 
the results of the present investigations and ends with an outlook for future 
studies. 



2 The nuclear energy density functional 

2.1 Framework and conjectures 

It is useful first to recall some basics of (non-relativistic) Density Functional 
Theory (DFT) [3,4,5,6] for fermionic many-body systems. The starting point 
of DFT is the Hohenberg-Kohn (HK) variational principle for the ground-state 
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energy of interacting fermions: 



E = min{F HK [p]} • 



(1) 



The exact ground-state energy is obtained from the universal functional Fhk 
depending only on the local density p(r). The Kohn-Sham DFT considers 
an auxiliary system of non-interacting quasi-particles and states that for any 
interacting system, there exists a local single-particle potential (the Kohn- 
Sham potential) such that the exact ground-state density of the interacting 
system equals the ground-state density of the auxiliary problem. The HK free 
energy is commonly decomposed into three separate terms: 



where is the kinetic energy of the non-interacting N-particle system, 
is a Hartree energy, and E xc denotes the exchange-correlation energy which, 
by definition, contains everything else. The practical usefulness of the Kohn- 
Sham scheme depends entirely on whether accurate approximations for E xc 
can be found. The local ground-state density is constructed using so-called 
auxiliary orbitals, 



which are unique functionals of the density p(r), i.e. the KS scheme defines 
a self-consistency problem. The set {4>ks} °f single-particle wave functions 
represents an "effective" basis because, in general, these functions do not have 
a directly observable physical interpretation 1 . The corresponding Kohn-Sham 
single-particle equations, which determine i^gg, are easy to solve numerically, 
even for large numbers of particles. The relativistic analogue of the Kohn- 
Sham scheme, used in the present work, starts with a set of Dirac spinors 
{ipn} which replaces the non-relativistic wave functions {ipxs}- 

The conjectures on which our approach to the nuclear energy density func- 
tional is based, with contact to low-energy QCD (see Sec. 1), can now be 
adapted as follows: 

(1) The large scalar and vector mean fields (with opposite signs) that have 
their origin in the in-medium changes of the chiral condensate and of the 
quark density, determine the Hartree energy functional -Eh[p]- 



except for the last occupied orbital (see Ref. [6]). 



F HK [p] = E kin [p] + E u [p] + E xc [p] , 




N 




(3) 
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(2) The chiral (pionic) fluctuations including one- and two-pion exchange 
with single and double virtual A(1232)-isobar excitations plus Pauli block- 
ing effects, determine the exchange-correlation energy functional E xc [p\. 

Note that we do not introduce explicit spin orbit terms. They emerge naturally 
from the Lorentz scalar and vector self-energies generated by the relativistic 
density functional, as discussed in the following section. Furthermore, contact 
terms from regularization-dependent pieces of two-pion exchange processes 
can be grouped together with the Hartree part E-r and need not be included 
in E xc . However, we prefer to keep them along with the pionic fluctuations for 
convenience. 



2.2 Density- dependent point coupling approach 



The density distribution and the energy of the nuclear ground state are ob- 
tained from self-consistent solutions of the relativistic generalizations of the 
linear single-nucleon Kohn-Sham equations. In order to derive those equations 
it is useful to construct a point-coupling model [18,19,20] with density depen- 
dent interaction terms, designed such as to reproduce the detailed density 
dependence of the nucleon self-energies resulting from E K [p] + E xc [p] . 

The spin-orbit force in nuclei is well known to be abnormally strong. An 
estimate based on the relation familiar from atoms, 

1 • s dU 

AHi s = 5 (4) 

2M 2 N r dr V ; 



using the average nuclear single particle potential U, would give a spin-orbit 
interaction too weak by about an order of magnitude. It is also known [17] 
that chiral one- and two-pion exchange dynamics alone cannot reproduce the 
observed 1 • s interaction. Additional mechanisms are required, such as the 
coherent action of strong scalar and vector potentials in relativistic mean-field 
models. In order to generate the large empirical spin-orbit splittings, we have 
two options. The first one is to stay within a non-relativistic framework (such 
as the Skyrme energy density functional) and simply fix the spin-orbit coupling 
by an appropriate contact term. The second option (the one we adopt here) 
is to preserve the distinction between scalar and vector nucleon self-energies 
in a relativistic description. 

In the effective field theory approach to the NN interaction [21] short range 
central and spin-orbit forces are represented by contact interactions with low- 
energy parameters fitted to NN phase shifts. It has been demonstrated in [22] 
how these low-energy parameters translate into an equivalent Skyrme repre- 
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sentation and - not surprisingly - account for the correct magnitude of the 
spin-orbit force as required already for the fit to triplet p-wave NN phase 
shifts. In our approach the role of these contact terms is taken over by the 
strong Lorentz scalar and vector fields which produce effectively the same 
spin-orbit coupling at the level of the nuclear mean field. 

A successful framework that meets these requirements for a two component 
system of protons and neutrons starts from a relativistic Lagrangian which 
includes isoscalar-scalar (S), isoscalar- vector (V), isovector-scalar (TS) and 
isovector-vector (TV) effective four-fermion interaction vertices with density- 
dependent coupling strengths. The density dependence represents many-body 
effects beyond mean field and two-body forces, including Pauli blocking effects. 
In Refs. [11,12] we have shown that a very good description of nuclear matter 
and finite nuclei from oxygen to calcium can be based on the following minimal 
Lagrangian density: 

(5) 



The four terms read: 



£fe» = $(i7/^-Mjv)V>, (6) 

-i G TS {p){^) • $7ty) - \ Gtv(p)(^t^) ■ (Hl^) , (7) 

*& } = ~ (8) 
C cm = eA*$^p->y l d> - \f, u F^ , (9) 

where ip is the Dirac field of the nucleon with its two isospin components 
(p and n). Vectors in isospin space are denoted by arrows. In addition to the 
free nucleon Lagrangian £f ree and the interaction terms contained in 4nt> when 
applied to finite nuclei, the model must include the coupling £ cm of the protons 
to the electromagnetic field A 11 with F^ v = d^Ay — d v A^ and a derivative 
(surface) term £ int . One could, of course, construct additional derivative terms 

(2) 

in £; nt , further generalized to include density dependent strength parameters. 
However, as we shall see, there appears to be no need in practical applications 
to go beyond the simplest ansatz (8) with a constant D s . In fact, present 
data on nuclear ground states constrain only this single isoscalar derivative 
term [23]. 

The variational principle 5£/5ip = applied to the Lagrangian (5) leads to 
the self-consistent single-nucleon Dirac equations, the relativistic analogue of 
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the (non-relativistic) Kohn-Sham equations. The nuclear dynamics produced 
by chiral (pionic) fluctuations in the medium is now encoded in the density 
dependence of the interaction vertices. 

In the framework of relativistic density functional theory [6,24,25,26], the 
density-dependent couplings are functions of the 4-current 

f = = faf , (10) 



where is the 4-velocity defined as (1 — v 2 ) 1 ' 2 (1, v). We work in the rest- 
frame of the nuclear system with v = 0. 

The couplings G;(p) {i = S, V,TS,TV) are decomposed as follows: 



G l (p)=Gf>+G ( ?\p) (fori = S,V) 
and d (p) = Gf ] (p) (for i — TS, TV) , (11) 

into density-independent parts G- ^ which arise from strong isoscalar scalar 
and vector background fields, and density-dependent parts G\ n \p) generated 
by (regularized) one- and two-pion exchange dynamics. As in our previous 
work, it is assumed that only pionic processes contribute to the isovector 
channels. 

The relativistic density functional describing the ground-state energy of the 
system can be re-written as a sum of four distinct terms: 

E [p] = E ircc [p] + E n [p] + E conl [p] +EJp], (12) 



with 



£frcc[p] = / d 3 r (M^PHl ■ V + M^lfo) , (13) 
E H [p] = \jdh {(0 o |4 o) (^) 2 |0o) + (0ol4 O) (^) 2 |0o)} , (14) 

E„[p} = \jd 3 r {(MG { s\p)(^) 2 \fo) + (<Po\G { v\p)(^) 2 \<Po) 

+ (0o|41(p)(^) 2 |0 o ) + (0ol4t(p)(^7M^) 2 |0o) 

- (0 o |4^[V(^)] 2 |0o)} , (15) 

E coul [p] = \jd?r (0 o |A^i^ 7 ^l0o) , (16) 

where \<f> ) denotes the nuclear ground state. Here E {ree is the energy of the 
free (relativistic) nucleons including their rest mass. E^i is a Hartree-type 
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contribution representing strong scalar and vector mean fields, later to be 
connected with the leading terms of the corresponding nucleon self-energies 
deduced from in-medium QCD sum rules. Furthermore, E n is the part of the 
energy generated by chiral 7riVA-dynamics, including a derivative (surface) 
term, with all pieces explicitly derived in [17]. 

Note that Eqs. (12-16) are arranged in close correspondence with a relativis- 
tic extension of the Hohenberg-Kohn energy functional (2), identifying the 
Hartree term E-& with Eq. (14) plus the Coulomb energy (16), and the ex- 
change correlation part with E n of Eq. (15). Minimization of the ground-state 
energy, represented in terms of a set of auxiliary Dirac spinors ipk , leads to the 
relativistic analogue of the Kohn-Sham equations. These single-nucleon Dirac 
equations are solved self-consistently in the "no-sea" approximation which 
omits the explicit contribution of the negative-energy solutions of the rel- 
ativistic equations to the densities and currents. This means that vacuum 
polarization effects are not taken into account explicitly. They are included 
in the adjustable parameters of the theory (for an exhaustive discussion see 
Refs. [24,27,28,29]). This no-sea approximation is well founded within the 
framework of an effective field theory. Nucleon-antinucleon fluctuations occur 
at high-energy (short-distance) scales, not resolved at the long wavelength 
characteristic of the nuclear Fermi momentum. Such fluctuations are implic- 
itly absorbed, together with other short-range dynamics, in contact terms of 
the effective Lagrangian. 

The expressions for the isoscalar and isovector four-currents and scalar densi- 
ties read: 



N 

3ii = (0o|^7^l0o) = Y^l^k , (17) 

k=l 

N 

i = (0o|^7 M ^l0o) = Y ^kl^Tlpk , (18) 

fc=l 

N 

Ps = (0o|^|0o) = Y^ipk, (19) 

k=l 

N 

Ps = (0o|^r^|0 o ) = Y^ripk, (20) 

fc=i 

where ipk are Dirac spinors and the sum runs over occupied positive-energy 
single-nucleon states. 

In this work we only consider systems with time-reversal symmetry in the 
ground-state, i.e. even-even nuclei. The space components of all currents vanish 
(j = 0), and because of charge conservation only the third component of isospin 
(r 3 = —1 for neutrons and r 3 = +1 for protons) contributes. The relevant 
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combinations of densities are: 



p=(0o|^7Vl0o)=P P + p n , (21) 

P3=(0o|^ 37 V|0o)=P P -p n , (22) 

Ps = (MM<M + (23) 

P53 = (0o|^|0o) = p^-pS- (24) 



The charge density is 

PcH = (0 O |^7°^^|0O) = P P ■ 

The single-nucleon Dirac equation is found by minimization with respect to 

• V + M N + 7 Sy + 7oT 3 £ T y 

+ 'j T, R + T, s + nT.Ts^k = , (25) 

with the self-energies: 





= 


E-ry 


— ^Ty 












■if 



+ Gl; ) (p)]p + eA i±^, (26) 

(p) Pa , (27) 
+ G% ) (p)]ps + D$ ) V 2 p a , (28) 
(P)P53, (29) 

^(p) 2 agW(p) , 

"dp - p + ^p— ps+ 

«9p Ps+ <9p P53 j ' (dUJ 

where A°(r) = ^/dV is the Coulomb potential. In addition to the 

usual contributions from the time components of the vector self-energies and 
the scalar potentials, we must also include the "rearrangement" terms, E/?, 
arising from the variation of the vertex functionals with respect to the nucleon 
fields in the density operator p. For a Lagrangian with density dependent cou- 
plings, the inclusion of the rearrangement self-energies is essential in order to 
guarantee energy-momentum conservation and thermodynamical consistency 
[30,31,32]. It is also required by the Hugenholtz-Van Hove theorem [33]. Using 
the single-nucleon Dirac equation and performing an integration by parts, the 
ground-state energy of a nucleus with A nucleons reads: 
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£o = X> -\j d 3 r { [Gf + Gf{p)] p% + 41(P) pf 3 + 

[4 0) + G , S; ) (p)]p 2 + 4t(p)P3 + 

dGf{ P ) dG? s ( P ) 2 

dp p + dp p3p + 
ep ch A° + DPp s V 2 p s } , (31) 

where denotes the single-nucleon energies. 



3 Linking the energy functional to the low-energy sector of QCD. 



3.1 Scalar and vector fields: guidance from in-medium QCD sum rules. 



The QCD ground state (or "vacuum") is characterized by strong conden- 
sates of quark-antiquark pairs and gluons, an entirely non-perturbative phe- 
nomenon. The quark condensate (qq), i.e. the ground state expectation value 
of the scalar quark density, plays a particularly important role as an order pa- 
rameter of spontaneously broken chiral symmetry. At a renormalization scale 
of about 1 GeV (with up and down quark masses m u +md ~ 12 MeV) the value 
of the chiral vacuum condensate [34] is (qq)o — —(240 MeV) 3 ~ —1.8 fm -3 . 
Hadrons, as well as nuclei, are excitations built on this condensed ground state. 
Changes of the condensate structure of the QCD vacuum in the presence of 
baryonic matter are a source of strong fields experienced by the nucleons in 
the bulk of nuclei. 

In- medium QCD sum rules relate the leading changes of the scalar quark 
condensate and of the quark density at finite baryon density, with the scalar 
and vector self-energies of a nucleon in the nuclear medium. To first order 
in the scalar and baryon densities, these self-energies can be expressed as 
follows [12,35,36,37,38,39]: 



4 o, = _<^ 



,( ) _ 4(m u + m d )M N 



where o~n = (N\m q qq\N) is the nucleon sigma term (~ 50 MeV), m n is the 
pion mass (138 MeV), and f n = 92.4 MeV is the pion decay constant. The 
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resulting E^ 1 and Ey are individually of the order of 300 - 400 MeV in 
magnitude. Their ratio 



a N PS 

4:(m u + m d ) p 



(34) 



is close to —1, suggesting a large cancellation of scalar and vector potentials in 
the single-nucleon Dirac equation, a feature characteristic of relativistic mean- 
field phenomenology. Of course the constraints implied by Eqs. (32)- (33), and 
by the ratio Eq. (34) are not very accurate on a quantitative level. Beyond the 
leading-order Ioffe formula [40] on which Eq. (32) is based, one finds correc- 
tions from condensates of higher dimension which are not well under control. 
Given in addition the uncertainties in the values of and m u + m^, the 
estimated error for the ratio Eg' l /Ey' ) ~ —1 is about 20%. It is important to 
note, however, that the explicitly calculated two-pion exchange fluctuations, 
used in this work, account for a large part of the effects from higher dimen- 
sional condensates (such as the four-quark condensate). Restricting E^y to 
the leading order terms (32-33) does therefore make sense. 

Comparing Eqs. (26) and (28) for the isoscalar vector and scalar potentials 
of the single-nucleon Dirac equations, with the Eqs. (32) and (33) for the 
condensate background self-energies, respectively, the following estimates hold 
for the couplings of the nucleon to the background fields (the Hartree terms 
in the energy functional): 

Gg» = , (35 ) 

TT'f TV 



and 



(0) _ 4(m u + m d )M N 
° V ~ </l • (36) 



Inserting the values of the nucleon and pion masses and the pion decay con- 
stant, 



4 0) ~-ll fm 2 



G{?~11 fm 2 



50 MeV 
A(m u + m d 
50 MeV 



(37) 
(38) 



which implies ~ —Gy 1 ~ —10.6 fm 2 for typical values crjv ~ 48 MeV and 
m u + rrid — 12 MeV (at a renormalization scale of about 1 GeV). This estimate 
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is useful for first orientation. In the actual applications to finite nuclei, G s v 
will have to be fine-tuned. We mention here in advance that, as a remark- 
able and non-trivial result of this fine-tuning, deviations from the "canonical" 
values (37-38) will turn out to be surprisingly small. 



3.2 The exchange- correlation term: In-medium chiral perturbation theory. 



The many-body effects represented by the exchange-correlation density func- 
tional are approximated by chiral 7riVA-dynamics, including Pauli blocking 
effects. Our framework is in-medium chiral perturbation theory with inclusion 
of A-isobars, computed to three-loop order in the energy density and described 
in detail in Refs. [13,14,15,16,17]. 

Regularization dependent short-range contributions from pion loops are en- 
coded in a few parameters representing counter terms, or equivalently, subtrac- 
tion constants in spectral representations of the respective momentum space 
amplitudes. At the level of the effective Lagrangian, they appear as contact in- 
teractions or derivatives thereof. The counter-term contributions to the energy 
per particle, E = E/A — M^, in symmetric nuclear matter are written: 

AE(k f )^ = 6 3 ^| + bS 4 + b 6 ^ . (39) 



At this stage the calculation for homogeneous nuclear matter excludes surface 
terms which are discussed in the next section. The constants foj are defined 
dimensionless, along with a convenient scale A to be chosen larger than the 
"small" scales (kf,m n and A) but smaller than the "chiral gap", 4irf n . Our 
choice is A = 2nf w ~ 0.58 GeV. 

For asymmetric nuclear matter, the energy per particle is written E as = E + 
S 2 S 2 with the asymmetry energy S 2 and 5 = (p n — ff)/p. The counter-term 
contributions to the asymmetry energy has an expansion analogous to (39), 

AS 2 (k f )^ = aS + aS 4 + , (40) 



with dimensionless constants aj. The parameters bi and a« subsume unresolved 
short-distance NN-dynamics in the nuclear medium. Note that no regulariza- 
tion is required at order kj where the corresponding loop integrals are all finite. 
There are no contact interactions that generate a kj term at the Hartree-Fock 
level; gradient terms start minimally as V 2 . At order the constants 65 and 
05 represent short- distance effects partly arising from momentum dependent 
interactions. In actual practice it turns out that these constants can be set 
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to zero [17], the primary momentum dependent forces being already well de- 
scribed by the finite parts of in-medium two-pion exchange pieces. Additional 
short-range three-body contributions and effects of higher loops which feed 
back into kj terms are parameterized by the constants b e and a 6 , respectively. 

In Ref. [17] the counter terms hi and ai have been adjusted to reproduce nuclear 
and neutron matter properties. The resulting equation of state of isospin- 
symmetric matter is in good agreement with recent microscopic calculations 
[41], though with a somewhat too high incompressibility (K = 304 MeV). The 
calculated real part of the nucleon single-particle potential, U(p,kf), is very 
close to the result of relativistic Dirac-Brueckner calculations (see Ref. [42]). 
Isospin properties of nuclear matter and the energy per particle in neutron 
matter are significantly improved by incorporating explicit 7riVA-dynamics, 
in comparison with earlier calculations (see Ref. [13]) which included only 
nucleons in two-pion exchange processes. 

While gross properties of infinite nuclear matter are useful for orientation, the 
large amount of nuclear observables studied in the present work provides a far 
more accurate data base that permits a better adjustment of the constants 
hi and a^. We refer to the more detailed discussion in Section 3.3 but point 
here already to the interesting result that the best fit values for 63 and 
turn out to be within only a few percent of those determined in the ChPT 
calculation of nuclear and neutron matter, while 65 = a 5 = can still be 
maintained. The only major difference is in the short-distance three-body 
term proportional to &6 f° r which the fit to the broad range of nuclear data 
requires stronger attraction. Throughout the procedure, the input values for 
the chiral pion-nucleon and ttNA couplings are strictly kept fixed by pion- 
nucleon scattering observables in vacuum. These couplings determine the finite 
parts of intermediate-range one- and two- pion exchange contributions to the 
energy density as predicted by in-medium ChPT. 

In the simplest DFT approach, the exchange-correlation energy for a finite 
system is determined in the local density approximation (LDA) from the 
exchange-correlation functional of the corresponding infinite homogeneous sys- 
tem, replacing the constant density p by the local density p(r) of the actual in- 
homogeneous system. In our case the exchange-correlation terms of the nuclear 
density functional are determined within LDA by equating the corresponding 
self-energies in the single- nucleon Dirac equation (25), with those arising from 
the in-medium chiral perturbation theory calculation of 7riVA-dynamics in ho- 
mogeneous isospin symmetric and asymmetric nuclear matter. Steps beyond 
the LDA will be taken by adding surface terms involving derivatives of the 
density. 

The density- dependent couplings are expressed as polynomials in frac- 
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tional powers of the baryon density: 



Gf \ P ) = a x + q 2P ^ + c aP i + c l4 p ... (i = S,V, TS, TV) . 



(41) 



The detailed derivation of the constants Cij is presented in the Appendix. 

The coefficient Dg of the derivative term in the equivalent point-coupling 
model (Eqs. (26) and (28)) can be determined from ChPT calculations for 
inhomogeneous nuclear matter. The density-matrix expansion method [43] is 
used to derive the in-medium insertion in the nucleon propagator, as shown 
in Refs. [14,17]. The isoscalar nuclear energy density emerging from chiral 
pion-nucleon dynamics has the form: 



The coefficients have been compared with the corresponding phenomenological 
parameters of various Skyrme type energy density functionals. In particular, 
the gradient term (Vp) 2 i<v(&;/) plays an important role in shaping the nu- 
clear surface [44]. While the applicability of the density- matrix expansion is 
questionable at very low densities, an important result of Refs. [14,17] is that 
around nuclear matter saturation density, the parameter-free ChPT predic- 
tion for Fy{kfo) is in very good agreement with the empirical values used in 
standard parameterizations of the Skyrme density functional [1]. We can ap- 
proximate Fy by a constant in the relevant region of densities. The following 
relationship between i*V and the derivative term of the point-coupling model 
(see Eqs. (8,15)) is then deduced by a straightforward comparison [45]: 



With the average value of F^(kf) in the region 0.1 fm 3 < p < 0.2 fm 3 [17], 
we estimate = — 0.7 fm 4 . 

The inclusion of derivative terms in the model Lagrangian and the determina- 
tion of its strength parameters from ChPT calculations in inhomogeneous mat- 
ter actually goes beyond the local density approximation. The term, Eq. (8), 
with the strength parameter given by Eq. (43), represents a second-order gra- 
dient correction to the LDA, i.e. the next-to-leading term in the gradient 
expansion of the exchange-correlation energy calculated by in-medium chiral 
perturbation theory. 



£(p,Vp)= P E(k f ) + (Vp) 2 F v (k f ) + ... . 



(42) 



- 2F V = D 



s 



(43) 
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Table 1 

Empirical nuclear matter properties (energy per particle, saturation density, incom- 
pressibility, Dirac effective mass, symmetry energy) employed in the adjustment of 
parameters. For each "empirical" quantity we include the range of values within 
which the quantity was allowed to vary during the fit. 





E/A (MeV) 


p sa t (fm 3 ) 


K (MeV) 


M*/M N 


S 2 (MeV) 


Fit 


-16. (±1.) 


0.15 (±0.01) 


250. (±25.) 


0.6 (±0.5) 


33. (±3.) 



Table 2 

Empirical properties of finite nuclei [57,58] employed in the fitting procedure: bind- 
ing energies (£*), charge radii (r c h) an d differences between r.m.s. neutron and pro- 
ton radii (r n — r p ). The relative weights used in fitting these quantities are 0.3%, 
0.2% and 10% respectively. 





E h (MeV) 


r c h (fm) 


r n ~ r p (fm) 


16Q 


-127.62 


2.73 




40 Ca 


-342.05 


3.49 




48 Ca 


-415.99 


3.48 




72 Ni 


-613.17 






90 Zr 


-783.89 


4.27 




116 Sn 


-988.68 


4.63 


0.12 


124 Sn 


-1049.96 


4.67 


0.19 


132 Sn 


-1102.92 






204p b 


-1607.52 


5.49 




208 pk 


-1636.45 


5.51 


0.20 


214p b 


-1663.30 


5.56 




210p o 


-1645.23 







3.3 Adjusting parameters 

The parameters of the point-coupling model are fixed simultaneously to prop- 
erties of nuclear matter (see Table 1), and to binding energies, charge radii and 
differences between neutron and proton radii of spherical nuclei (see Table 2). 
The inclusion of neutron-rich nuclei, in particular, provides strict constraints 
for the parameters in the isovector channel. In both tables the empirical val- 
ues and data are listed together with the relative errors used in the fitting 
procedure. 
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Corrections have been applied for center-of-mass motion and pairing. From 
the solution of the self-consistent equations we subtract the microscopically 
calculated center-of-mass correction from the total binding energy, following 
[46]: 

(P 2 ) 

AE cm = - V , (44) 
2AM N ' v ; 



where P c 2 m is the total squared momentum of a nucleus with A nucleons. For 
the open shell nuclei in Table 2, pairing correlations are treated in the BCS 
approximation with empirical gaps [47]. 

We have adjusted a minimal set of parameters starting from the estimates Eqs. 
(37) and (38) for the couplings of the condensate background fields (Hartree 
term), and the counter term constants b 3 , a 3 , b 5 , a 5 and b 6 , a 6 guided by 
Ref. [17] for the self-energies arising from chiral 7riVA-dynamics (exchange- 
correlation term). In particular we set 65 = 05 = as in Ref. [17]. The total 
number of adjustable parameters is seven. However, five of those will turn out 
to be remarkably close to preceding expectations or predictions. Only the two 
constants associated with three-body correlations are free of links to other 
constraints. 

The resulting optimal parameter set (FKVW) is shown in Table 3 in compar- 
ison with estimates and predictions: a) for the Hartree term from in-medium 
QCD sum rules; b) for the exchange-correlation functional from the in- medium 
ChPT calculations of 7riVA-dynamics in nuclear matter. 

Evidently the best fit parameters are remarkably close to the anticipated QCD 
sum rule and ChPT values, with the exception of the /cj terms (b^ and a^) 
for which the fit to nuclear data systematically requires an attractive shift as 
compared to the ChPT calculation [17]. 

The QCD sum rule predictions for the condensate background scalar and vec- 
tor self-energies have an accuracy of about 20 % as mentioned previously. Sur- 
prisingly, however, the leading-order sum rule estimates, Eqs. (37) and (38), 
turn out to be realized remarkably well when the parameters of the scalar 
and vector Hartree fields are varied freely and optimized in confrontation 
with a large number of high-precision nuclear data. There appears to be an 
a posteriori consistency with the implicit assumption made in the in-medium 
ChPT calculations based on 7riVA-dynamics: namely that the individually 
large scalar and vector nucleon self-energies Eg and E^/ 1 cancel in their con- 
tribution to the energy density around saturation. Of course, the 63 counter 
term at order from regularized two-pion exchange can be translated into 
an equivalent Hartree term as well, such that the total balance of scalar and 
vector self-energy contributions linear in density is attractive and reminiscent 
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Table 3 

Best-fit parameters of the density-dependent point coupling model (Fit FKVW). 
Upper panel: coupling strengths of background (Hartree) terms in comparison with 
estimates from in- medium QCD sum rules (37, 38) (QCDSR). Lower panel: con- 
stants of counter terms from short-distance regularization (reference scale A = 

(it) 

2irfi T ~ 0.58 GeV) and surface (derivative) term D s in comparison with re- 
sults from in-medium ChPT at three-loop order including A(1232) [17]. The choice 
^5 = 05 = from Ref. [17] has been kept, and these parameters are not included in 
the fitting procedure. 



Background (Hartree) terms 






Fit FKVW 




QCDSR 






- 11.5 




-10.6 






11.0 




10.6 


Counter terms / Pionic fluctuations 






Fit FKVW 




ChPT [17] 






-2.93 




-3.05 






2.20 




2.16 


K 
























&6 




-5.68 




-2.83 


«6 




-0.13 




2.83 


[fm 4 ] 




-0.76 




-0.7 


*: inputs, see Ref. [17] 



of QHD and the Walecka model [27]. The repulsive term of order kj (or p 4 ^ 3 ), 
on the other hand, results model-independently from Pauli blocking effects on 
chiral two-pion exchange within in-medium ChPT and has no counterpart in 
relativistic mean field models. 

It is instructive to examine the hierarchy of terms in the /^-expansion of the 
nucleon isoscalar single particle potential, 

k 3 k 4 k 5 k 6 k 7 
U(k f ) = ^ + 9^+9^ + 9^ + 0^) , (45) 

given in detail in the Appendix. The reference scale is chosen again as A = 
2nf n ~ 0.58 GeV. The resulting coupling strengths g t are listed in Table 
4. They are all of "natural" size when looked at from the point of view of 
an effective field theory. Background fields can be ignored here because their 
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Table 4 

Comparison between fine-tuned coefficients gi (Fit FKVW) in the expansion (45) 
of the isoscalar single-particle potential, and ChPT values of g\ deduced from cal- 
culations performed in Ref. [17]. The constants 6, used in this compilation are those 
given in Table 3 (see Appendix for further details). 





Fit FKVW 


ChPT [17] 


93 


-2.84 


-3.10 


94 


2.67 


2.67 


95 


2.17 


2.17 


96 


-3.35 


5.21 



contributions to U cancel almost completely, i.e. = Gg ps + Gy'p ~ 0. 

An alternative way of displaying the convergence properties of the kf ex- 
pansion is by translating U of Eq. (45) into the following form of a (local) 
potential: 

U{I) = U °W) +U '{W)) +Ui {m) +Ue {W)) ■ (46) 



where p(0) is the central (r = 0) value of the density distribution p(r). As 
a typical example (for p ~ p(0) ~ p sat ) one finds a sequence of coefficients 
U n as shown in Fig. 1. Note once again that the leading attractive piece U3 
(proportional to the density p) does not receive contributions from the individ- 
ually much stronger background self-energies £y ~ — Eg ~ 0.35 GeV which 
cancel in the energy per particle but coherently build up the large spin-orbit 
coupling in finite nuclei. The terms of order p 4 / 3 and p 5 / 3 (with coefficients C/4 
and C/5) are repulsive and successively smaller, while the "three-body" term 
Uq proportional to p 2 is attractive and even smaller in magnitude. 

It might appear that the attractive p 2 term causes nuclear matter to collapse 
at very high density. We have checked that stability persists up to at least 
three times the density of normal nuclear matter. At even higher densities, 
expansions such as (45) must be carried to higher orders in kf. One should 
recall, however, that the limit of applicability for the present chiral approach 
is imposed by the condition that kf must remain sufficiently small compared 
to the reference scale A. 

Using the input parameters specified earlier in this section, let us have a 
look at a first series of results achieved in comparison with data. In Fig. 2 
we plot the relative deviations between calculated and experimental values, 
{{O l th — Ol )/O l e ) in percent, of binding energies and charge radii for all 
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nuclei that have been included in the fit (see Tab. 2). The RMF+BCS results 
obtained with the FKVW parameter set are compared to those of two standard 
phenomenological relativistic mean-field models using a) the non-linear meson 
exchange interaction NL3 [48], and b) the density-dependent meson-exchange 
interaction DD-ME1 [49]. The binding energies and charge radii calculated in 
the present work are comparable or superior to those produced with two of 
the most accurate phenomenological relativistic mean-field interactions (the 
number of free parameters is seven for NL3, and eight for DD-ME1). 



4 Ground-state properties of spherical and deformed nuclei 

In Ref. [12] a first version of our relativistic point-coupling model (at that stage 
without explicit inclusion of A(1232) degrees of freedom) has been tested in 
an analysis of the equations of state for symmetric and asymmetric nuclear 
matter, and of bulk and single-nucleon properties of light and medium heavy 
N ~ Z nuclei. The detailed isospin dependence of the effective interaction 
was, however, less than perfect. For heavier N > Z nuclei such as 208 Pb 
the calculated binding energies were deviating considerably (by more than 5 
%) from the experimental values. In the present version of the model, both 
the isoscalar and the isovector channels of the exchange-correlation functional 
incorporate the additional contribution of two-pion exchange with single and 
double virtual A(1232)-isobar excitations. Much better isospin properties of 
nuclear matter have thus been obtained. One can therefore expect improved 
results also for ground-state properties of N ^ Z nuclei calculated with the 
advanced nuclear energy density functional derived the present work. This is 
illustrated in Fig. 3 where for the set of spherical nuclei that were used in the fit 
(see Tab. 2), we compare the deviations between calculated and experimental 
binding energies and charge radii, respectively, for the FKVW interaction with 
explicit inclusion of A(1232) excitations, and for the previous version of the 
interaction [12] which did not include the A(1232) degree of freedom. The 
latter displays a systematic underbinding with the increase of the neutron - 
proton asymmetry in heavier nuclei, and the deviations exceed 5% already for 
Sn isotopes. This is corrected with the additional contribution of the A(1232)- 
isobar, and the new effective interaction reproduces with high accuracy the 
empirical masses. We also notice an improvement for the calculated charge 
radii (lower panel of Fig. 3), although not as dramatic as in the case of binding 
energies. 

In this section the improved effective FKVW interaction will be tested in self- 
consistent calculations of ground-state observables for spherical and deformed 
medium-heavy and heavy nuclei. The calculations, including open-shell nuclei, 
are performed in the framework of the relativistic Hartree-Bogoliubov (RHB) 
model, a relativistic extension of the conventional Hartree-Fock-Bogoliubov 
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method, that provides a basis for a consistent microscopic description of 
ground-state properties of medium-heavy and heavy nuclei, low-energy ex- 
cited states, small-amplitude vibrations, and reliable extrapolations toward 
the drip lines [2]. In the particle-hole channel we employ the new microscopic 
FKVW interaction. In comparison we also present calculations, for spheri- 
cal nuclei, using the DD-ME1 interaction. The DD-ME1 model Lagrangian is 
characterized by a phenomenological density dependence of the a, uj and p 
boson-nucleon vertex functions, adjusted to properties of nuclear matter and 
finite nuclei. The DD-ME1 effective interaction has recently been employed 
in a number of studies of ground-state properties of spherical and deformed 
nuclei [49,50], and of multipole giant resonances in spherical nuclei [51,52,53]. 
When compared to results obtained with standard nonlinear relativistic mean- 
field effective forces, the DD-ME1 interaction gives an improved description of 
asymmetric nuclear matter and of ground-state properties of N ^ Z nuclei. In 
the relativistic random phase approximation (RRPA), the DD-ME1 effective 
interaction reproduces the experimental excitation energies of multipole giant 
resonances in spherical nuclei. 

Pairing effects in nuclei are restricted to a narrow window of a few MeV 
around the Fermi level. Their scale is well separated from the scale of binding 
energies which are in the range of several hundred to thousand MeV. There 
is no empirical evidence for any relativistic effect in the nuclear pairing field. 
Therefore pairing can be treated as a non- relativistic phenomenon [54]. In 
most applications of the RHB model the pairing part of the well known and 
successful Gogny force [55] has been used in the particle-particle channel. We 
will conveniently proceed in the same way, noting at the same time that the 
chiral N-N potential produces very similar pairing matrix elements [56]. 



4-1 Spherical nuclei 

In order to test the detailed isospin dependence we investigate now the system- 
atics of isotopic chains. Deviations (in percent) of calculated binding energies 
from their experimental values [57] are shown in Fig. 4 for a series of even-A Sn 
isotopes. In addition to the FKVW interaction, computations are performed 
with phenomenological density-dependent interaction DD-ME1. The Gogny 
D1S force is used for the pairing interaction. Both with FKVW and DD-ME1, 
very good results are found over the entire major shell 50 < N < 82. For 
the new microscopic FKVW interaction in particular, the absolute deviations 
of the calculated masses from data do not exceed 0.1 %. In the lower panel 
of Fig. 4 we display the calculated charge radii of Sn isotopes in comparison 
with the experimental values [58]. For both interactions the theoretical values 
are in excellent agreement with data. 
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The isotopic dependence of the deviations (in percent) between the calculated 
binding energies and the experimental values for even-A Pb nuclei [57], is 
plotted in the upper panel of Fig. 5. It is interesting to note that, although 
DD-ME1 and FKVW represent different physical models, they display a sim- 
ilar mass dependence of the calculated binding energies for the Pb isotopic 
chain. On a quantitative level the FKVW interaction produces better results, 
with the absolute deviations of the calculated masses below 0.1 % for A > 190. 
In lighter Pb isotopes one expects that the observed shape coexistence phe- 
nomena will have a pronounced effect on the measured masses. 

Because of the intrinsic isospin independence of the effective single-nucleon 
spin-orbit potential, relativistic mean-field models naturally reproduce the 
anomalous charge isotope shifts [60]. The well known example of the anoma- 
lous kink in the charge isotope shifts of Pb isotopes is illustrated in the lower 
panel of Fig. 5. The results of RHB calculations with the DD-ME1 and FKVW 
effective interactions are shown in comparison with experimental values [61]. 
Both interactions reproduce in detail the A-dependence of the isotope shifts 
and the kink at 20S Pb. 

The differences between the r.m.s. radii of neutron and proton ground-state 
distributions of Sn and Pb nuclei are shown in Fig. 6. Results obtained with 
the FKVW interaction are compared with data [62] for the Sn isotopes, and 
with the empirical value of r n — r p in 208 Pb (0.20 ± 0.04 fm from proton 
scattering data [63], and 0.19 ± 0.09 fm from the excitation of the isovector 
giant dipole resonance by a-scattering [64]). 

The determination of neutron density distributions provides not only basic 
nuclear structure information, but it also places important additional con- 
straints on the isovector channel of the effective interactions used in nuclear 
models. In a recent analysis of neutron radii in non-relativistic and covari- 
ant mean-field models [65], the linear correlation between the neutron skin 
and the asymmetry energy was studied. In particular, it has been shown that 
there is a very strong linear correlation between the neutron skin thickness in 
208 Pb and the parameters that determine the asymmetry energy. The excel- 
lent agreement between the calculated r n — r p and the available data confirms 
that the isovector channel of the microscopic effective interaction FKVW is 
correctly represented. Note that this is achieved without ad-hoc introduction 
of a phenomenological p-meson exchange interaction. 

The FKVW model is based on the framework of density functional theory 
and, therefore, a natural test is the comparison of the calculated density dis- 
tributions with available data. In Fig. 7 we display several calculated charge 
form factors in comparison with data [66]. In addition to 48 Ca, 90 Zr, and 
208 Pb, for which the charge radii i.e. the first minima of the charge form fac- 
tors have been used in the fit of the parameter set FKVW (see Table 2), we 
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also plot the charge form factors of 92 Mo, 94 Zr, and 144 Sm. The theoretical 
charge density distributions are obtained by folding the point proton density 
distribution with the proton charge distribution of exponential form ~ e~ Ar 
(reflecting a dipole form factor), with the empirical proton charge radius. The 
corresponding form factors are plotted as functions of the momentum transfer 

q 



where p c h{o) is the Fourier transform of the spherical charge density and the 
second term is a center-of-mass correction. Higher-order effects in q 2 /P^ m are 
negligible. The excellent agreement between the calculated and experimental 
charge form factors for momenta \q\ < 2.5 fm -1 demonstrates that the FKVW 
interaction reproduces not only the moments of the distributions, but also the 
detailed charge density profiles. 

One of the principal advantages of using the relativistic framework lies in 
the fact that the effective single-nucleon spin-orbit potential arises naturally 
from the Dirac equation. The single-nucleon potential does not introduce any 
adjustable parameter for the spin-orbit interaction. In the FKVW model, in 
particular, the large effective spin-orbit potential in finite nuclei is generated 
by the strong scalar and vector condensate background fields of about equal 
magnitude and opposite sign, induced by changes of the QCD vacuum in 
the presence of baryonic matter [12]. In Fig. 8 we plot the deviations (in 
percent) between the calculated and experimental values of the energy spac- 
ings between spin-orbit partner-states in a series of doubly closed-shell nuclei. 
Even though the Kohn-Sham orbitals cannot be directly identified with the 
single-nucleon states in the interacting system, the calculated energy spacings 
between spin-orbit partners, especially the ones close to the Fermi surface, 
compare well with the experimental values. The experimental data are from 
Ref. [67], and the theoretical spin-orbit splittings have been calculated with 
the FKVW and DD-ME1 interactions. For the phenomenological DD-ME1 in- 
teraction the large scalar and vector nucleon self-energies which generate the 
spin-orbit potential, arise from the exchange of "sigma" and "omega" bosons 
with adjustable strength parameters. We notice that, even though the values 
calculated with DD-ME1 are already in very good agreement with experi- 
mental data, a further improvement is obtained with the FKVW interaction. 
This remarkable agreement indicates that the initial estimates Eqs. (37) and 
(38) for the condensate background couplings have perhaps been more realis- 
tic than anticipated, considering the uncertainties of lowest-order in-medium 
QCD sum rules. 




(47) 
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4-2 Deformed nuclei 



Deformed nuclei with N > Z present further important tests for nuclear struc- 
ture models. Ground-state properties, in particular, are sensitive to the isovec- 
tor channel of the effective interaction, to the spin-orbit term of the effective 
single-nucleon potentials and to the effective mass. 

In this section we test our advanced model in the region 60 < Z < 80. We 
compare predictions of the RHB calculations for the total binding energies, 
charge radii and ground-states quadrupole deformations of nine even-Z iso- 
topic chains with available data. The FKVW effective interaction is used in 
the particle-hole channel, and pairing correlations are described by the pairing 
part of the finite range Gogny D1S interaction [55]. 

The deviations (in percent) of the calculated binding energies from the exper- 
imental data [57] for Nd, Sm, Gd, Dy, Er, Yb, Hf, W, and Os isotopes are 
plotted in Fig. 9. Good agreement is found over the entire region of deformed 
nuclei. The maximum deviation of the calculated binding energies from data 
is below 0.5% for all isotopes. 

In Fig. 10 we compare the calculated charge radii with data from Ref . [58] . The 
charge density is constructed by folding the calculated (point) proton density 
distribution with the empirical proton-charge distribution. The calculation of 
the charge rms radius from the full charge form factor is quite involved for 
deformed nuclei [1,59]. We employ the frequently used simplified expression 



where r 2 is the nuclear mean-square radius of the (point) proton density dis- 
tribution. The calculated charge radii reproduce in detail the experimental 
isotopic trends. 

The ground-state quadrupole deformation parameter /?2, proportional to the 
expectation value of the quadrupole operator ((f> \3z 2 — r 2 \(f> ), is calculated ac- 
cording to the prescription of Ref. [68]. The theoretical values of the quadrupole 
deformation parameter are displayed in Fig. 11, in comparison with the em- 
pirical data extracted from B(E2) transitions [69]. We notice that the RHB 
results reproduce not only the global trend of the data but also the saturation 
of quadrupole deformations for heavier isotopes. 

The structure of nuclei far from stability provides a particularly important 
testing ground for global effective interactions. An interesting question, there- 
fore, is how far from stability can we extrapolate the predictions of the FKVW 
interaction. While an extensive analysis of this problem is beyond the scope 




(48) 
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of the present work, we consider in detail the example of neutron-deficient Pb 
isotopes which exhibit an interesting variety of coexisting shapes. They are 
characterized by the competition of low-lying prolate and oblate minima. Re- 
cent data [70] on energy spectra and charge radii indicate that, although the 
ground states of neutron-deficient Pb nuclei are spherical, both oblate (two- 
particle - two-hole proton excitations across the Z = 82 major shell) and 
prolate (four-particle - four-hole proton excitations) low-lying minima can be 
observed for N < 110 at almost identical excitation energies. 

A number of theoretical analyses of shell quenching and shape coexistence 
phenomena in neutron-deficient H g, Pb and Po nuclei has been reported 
in recent years. The quantitative description of coexisting spherical, oblate 
and prolate minima necessitates restoration of broken symmetries and the 
explicit treatment of quadrupole fluctuations. In particular, excellent results 
for shape coexistence in neutron-deficient Pb isotopes have recently been ob- 
tained by performing configuration mixing of angular momentum projected 
self-consistent mean-field states, calculated with the finite range and density 
dependent Gogny interaction [71] , and with the Skyrme interaction SLy6 sup- 
plemented by a density-dependent zero-range pairing force [72]. 

The model employed in the present work does not include angular momentum 
projection, nor can it account for configuration mixing effects. Nevertheless, 
we can compare the potential energy curves, calculated as functions of the ax- 
ial quadrupole deformation, with the corresponding mean-field energy curves 
calculated with the Gogny [71] and Skyrme [72] interactions. 

In Fig. 12 we display the calculated binding energy curves of even-A 182 - 196 pfr 
isotopes as functions of the quadrupole deformation. The curves correspond 
to axially deformed RHB model solutions with constrained quadrupole de- 
formation. The effective interaction is FKVW + Gogny D1S (pairing). The 
general agreement with the mean-field potential energy curves of Refs.[71] 
and [72] is satisfactory and, in particular, the present calculation predicts the 
coexistence of oblate and prolate minima at almost identical excitation ener- 
gies in 186 > 188 P5 ; in agreement with data. The relative excitation energies of 
coexisting minima based on different intruder configurations are determined 
by the spherical magic and semi-magic energy gaps. For A > 190 the mini- 
mum on the prolate side vanishes, in accordance with the nonrelativistic cal- 
culations performed with the Gogny and Skyrme-SLy6 interactions. In these 
nuclei, however, we also calculate a tiny shift of the ground-state minimum 
toward prolate deformation, a result not predicted by the two non-relativistic 
interactions, and not corroborated by available data. Although the quanti- 
tative analysis of the structure of neutron-deficient Pb nuclei must include 
additional correlations related to the restoration of broken symmetries and 
quadrupole fluctuations, it is indeed an encouraging result that the FKVW 
interaction reproduces the isotopic trend of the coexistence of shapes in heavy 
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nuclei far from stability. This reinforces the confidence in the detailed behav- 
ior of the isospin-dependent interaction derived from chiral two-pion exchange 
dynamics, improved by explicit A(1232) components. 



5 Summary, conclusions and further comments 

We have derived a relativistic nuclear energy density functional with connec- 
tions to two closely linked features of QCD in the low-energy limit: 

a) in-medium changes of vacuum condensates; 

b) spontaneous chiral symmetry breaking. 

The leading changes of the chiral (quark) condensate and quark density in 
the presence of baryonic matter are sources of strong (attractive) scalar and 
(repulsive) vector fields experienced by nucleons in the nucleus. These fields 
produce Hartree potentials of about 0.35 GeV in magnitude at nuclear matter 
density, in accordance with QCD sum rules. While these scalar and vector 
potentials cancel approximately in their contribution to the energy, they are 
at the origin of the large spin-orbit splitting in nuclei. 

The spontaneously broken chiral symmetry in QCD introduces pions as Gold- 
stone bosons with well-defined (derivative) couplings to baryons plus symme- 
try breaking corrections. In the present approach the exchange-correlation part 
of the energy density functional is deduced from the long- and intermediate- 
range interactions generated by one- and two-pion exchange processes. They 
have been computed using in-medium chiral perturbation theory with explicit 
inclusion of A(1232) degrees of freedom which turn out to be important. Regu- 
larization dependent contributions to the energy density, calculated at three- 
loop level, are absorbed in contact interactions with constants representing 
unresolved short-distance dynamics. 

This framework is translated into a point-coupling model with density-dependent 
interaction vertices. This is done for the practical purpose of deriving and 
solving self-consistent Dirac equations (relativistic analogues of Kohn-Sham 
equations) in order to determine the nucleon densities which enter the energy 
functional. We demonstrate that this scheme works extremely well when con- 
fronted with a large number of high-precision nuclear data over a broad range 
of spherical and deformed nuclei. In fact the quantitative accuracy of the cal- 
culated binding energies and radii is such that deviations from empirical data 
are generally less than 0.5% throughout the nuclear chart and usually much 
less (0.1%) for medium-heavy nuclei. It has to be noted, however, that these 
results still do not reach the accuracy of the best phenomenological mass ta- 
bles (recent trends in the determination of nuclear masses have been reviewed 
in Ref. [73]). Modern Skyrme-based microscopic mass formulas, with a total 
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of ~ 20 empirical parameters, fit the measured masses of more than two thou- 
sand nuclei with an rms error of less than 700 keV. In addition, they also 
produce excellent results for charge radii, with an rms deviation of ~ 0.025 fm 
for the absolute charge radii and charge isotope shifts of more than 500 nuclei. 
The nuclear energy density functional developed in this work is not, of course, 
a mass formula and, with only a small set of adjustable parameters, is not 
designed to compete with such a level of accuracy. Nevertheless, it provides 
an excellent microscopic framework on which a more fundamental approach 
to the calculation of nuclear masses can be based. 

A special focus in this work is on systematic studies along isotopic chains. 
Binding energies, as well as proton and neutron radii of such systems, with 
fixed number of protons but varying number of neutrons, are an excellent 
testing ground for the detailed isospin properties of the effective interaction. 
The quality of the results points to the fact that chiral pion dynamics, and 
especially the Van der Waals-like two-pion exchange forces with inclusion of 
the A(1232), is capable of generating the proper isospin dependence required 
by observables in isotopic sequences. 

The construction of the density functional involves an expansion of nucleon 
self-energies in powers of the Fermi momentum up to and including terms of 
order k®, or equivalently, 0(p 2 ) in the proton and neutron densities. Up to 
this order the present model has seven parameters, four of which are related 
to contact (counter) terms that appear in the chiral perturbation theory treat- 
ment of nuclear matter. One parameter fixes a surface (derivative) term and 
two more represent the strengths of scalar and vector Hartree fields. 

In the "best fit" set which reproduces a large amount of data on nuclear 
ground state properties, five of those seven parameters turn out to be sur- 
prisingly close to estimates and predictions from in-medium QCD sum rules 
and ChPT calculations for nuclear matter. In particular, the strong scalar and 
vector (Hartree) fields required to reproduce spin-orbit splittings are remark- 
ably compatible with the leading-order QCD sum rule estimates. The fitted 
coefficient 63 multiplying the O(k^) term in the chiral expansion of the nucleon 
self-energy, is within less than 5 % of the one from the ChPT calculation for 
nuclear matter. For the corresponding coefficient 03 in the expansion of the 
isospin dependent nucleon self-energy, the deviation (between fit to the nu- 
clear data base and ChPT evaluation of the symmetry energy) is even less 
than 2 %. The surface term required to reproduce the detailed systematics of 
nuclear radii is within less than 10 % of the ChPT prediction for slightly in- 
homogeneous nuclear matter. This prediction is actually parameter-free since 
the corresponding loop integrals involve only finite, regularization indepen- 
dent pieces. The exceptions that deviate from this overall pattern are the two 
parameters related to the three-body contact interaction terms (of order k®, 
or p 2 , in the self-energies) for which the nuclear data require more attraction 
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than anticipated by the in-medium ChPT estimate. On absolute scales, how- 
ever, these three-body terms are relatively small within the hierarchy of the 
kf expansion: an encouraging result although it cannot strictly be taken as 
proof of convergence. 

In this work we have only considered the low-energy QCD constraints on the 
effective interaction in the particle-hole channel, whereas the pairing part of 
the Gogny D1S force has been used for the effective interaction in the particle- 
particle channel. One could, of course, extend the chiral EFT approach to the 
pairing channel, and this would introduce additional contact terms. The first 
step in this direction has been taken in Ref. [56], where an improved version of 
the chiral nucleon-nucleon potential at next-to-next-to-leading order has been 
used to calculate the pairing gap in isospin-symmetric nuclear matter. The 
pairing potential consists of long-range one- and two-pion exchange terms, and 
two NN-contact couplings. It has been shown that the inclusion of the two- 
pion exchange at the next-to-next-to-leading order reduces substantially the 
cut-off dependence of the ^So pairing gap, indicating reasonable convergence of 
the small momentum expansion. The results are very close to those obtained 
with the phenomenological Gogny D1S force. 

The exchange-correlation energy functional in nuclear matter, calculated by in- 
medium chiral perturbation theory, has been used in Kohn-Sham calculations 
of finite nuclei by employing a second-order gradient correction to the local 
density approximation. Even though the quality of the results is on the level 
of the best phenomenological (non-relativistic and relativistic) self-consistent 
mean-field models, obviously the goal is to further improve the accuracy of the 
calculated nuclear ground-state energies and density distributions all over the 
periodic table. In the next step this will require the development of an accu- 
rate generalized gradient approximation for the nuclear exchange-correlation 
energy. 

An important issue will also be the relation with the universal low-momentum 
interaction Vi ow -k [74], deduced from phase-shift equivalent nucleon-nucleon 
interactions. An apparent close correspondence between perturbative in-medium 
chiral dynamics and Vi ow -k in the low-density limit has already been pointed 
out recently in Ref. [17]. This is a non-trivial observation, given that the input 
to the chiral approach [17] has been fixed to nuclear matter properties and not 
to NN phase shifts. The applicability of a perturbative expansion scheme for 
nuclear matter with low-momentum interactions, which is the basic hypoth- 
esis behind our approach, is discussed in Ref. [75]. Further studies towards 
foundations along these lines are presently being pursued. 

This work has demonstrated that chiral effective field theory provides a con- 
sistent microscopic framework in which both the isoscalar and isovector chan- 
nels of a universal nuclear energy density functional can be formulated. The 
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present approach to nuclear DFT establishes a fundamental link between low- 
energy QCD and ground-state properties of finite nuclei. From a practical 
point of view, a fully microscopic basis of effective nuclear interactions is es- 
pecially important for studies of nuclear structure in regions far from the valley 
of /3-stability, where extrapolations of phenomenological (non-relativistic and 
relativistic) models lack predictive power. 
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7 Appendix 



Here we present details of the nucleon self-energies from in-medium chiral 
perturbation theory and their translation into a point- coup ling model with 
density-dependent vertices. In [13,14,17] the momentum and density-dependent 
single particle potential of nucleons in asymmetric nuclear matter: 

UfakA-UfakffoS + OiP) with 5=^^, (49) 



has been calculated within two-loop in-medium chiral perturbation theory 
(ChPT). 

From the sums 



U(p = k f ,k f ) = ^ T (k f ,p) + Z™ >T (k f ,p) (50) 
-Urip = k ft k f )6 = Z%f T (k f , p) + Z™ PT (k f , p) 

of the scalar and vector nucleon self-energies in the isoscalar and isovector 
channels can be extracted at the Fermi surface, p = kf. It has been shown 
that the differences between scalar and vector parts of the ChPT self-energies 
are very small (suppressed by factors M^ 2 ) and linearly proportional to p. 
Their effects can therefore be absorbed as small corrections to the condensate 
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fields. This motivates the following ansatz: 



^ChPT 



(k f ,p) = W(k f ,k f ) 



^{k f ,p) = \U{k f ,k f ) 

( 51 ) 

?$F{k f ,p)=-\U I {k f ,k f )5 
Tffi £ {kt iP ) = -\U I {k h kt)6 



in which the single-nucleon self-energy Eq. (50) is equally divided between the 
vector Eq. (26) and the scalar Eq. (28) components. For the ChPT isoscalar 
and isovector self-energies [17] a polynomial fit up to order k 6 f is performed: 



U(k f ;b 3 ,b 5 ,b 6 ) 



k 3 k 4 



+ h + 36 6 ] ^ 



0{kl) 



C 5 + 



1 *J 
A 4 



(52) 



and 



Ui(kf, a 3 , a 5 , a 6 , b 5 ) — [d 3 + 2a 3 ] -± 



k 4 

+ d4 Ts + 

k 6 

+ [d 6 + 2a 6 ]^ 



d 5 + 2a 5 - —6 5 

y 



A 4 



(53) 



where &j and a« are the dimensionless short-distance regularization constants 
appearing in the counter terms (39) and (40). They subsume all unresolved 
short-range two- and three-body dynamics, plus possible higher order (four- 
loop etc.) effects not evaluated but contributing already at order kf. The 
constants q and d% are derived directly from the in-medium ChPT calculations 
of the finite (regularization-independent) parts of one- and two-pion-exchange 
contributions to the energy, with inclusion of A(1232) intermediate states. 
They are taken over unchanged from Ref. [17]. Note that the coupling strengths 
gi in the isoscalar single-particle potential (45) are identified as 



g 3 = c 3 + 2b 3 . 

94 = C4 , 

#5 = c 5 + -65 , 
5 



96 = c e 



36 fi . 



(54) 
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The counter term constants bi and a, are listed in Table 3. The constants q 
and di are collected in Table 5 at the reference scale A = 2nf n ~ 0.58 GeV. 

In order to determine the density- dependent couplings of the exchange corre- 
lation pieces generated by the point-coupling model, the ChPT self-energies 
are first re-expressed in terms of the baryon density p = 2fcj/37r 2 = pP + p n 
and the isovector density p 3 = p p — p n : 



1 2 



E s hPT (^/> P) = ( C sl + c s2p^ + c s3 p f + c s4 p) p , (55) 
^v hPT ( k fi P) = (°vl + Cv2p T > + c v3 pt + c vi p) p , (56) 
S T5 PT ( /c / ) P) = (°tsl + Cts2p~ z + CtsZp* + C tsi p) p 3 , (57) 

E^ PT (/c / , p) = (c to i + c tv2 p3 + c to3 p3 + c tvi p) p 3 , (58) 

Next these self-energies are identified with the corresponding contributions to 
the potentials in the point coupling single-nucleon Dirac equation (26)-(30): 



T^(k,,p)=G^ ) {p)p 8 , (59) 
E^ / ,p)^«(p)p + i{^/ + ^pl} , (60) 

T$%{k ft p)=G£l(p) Pia , (61) 
ti%(k f ,p)=G$(p)p 3t (62) 

The resulting expressions for the density-dependent couplings of the pionic 
fluctuation terms read 2 



( \ 12 

G s (p) = °si + c s2p* + c s3 p 3 + c s4 p , (63) 

( \ 12 

G ( y> (p) = c vl + c v2 p 3 + c v3 p3 + c v4 p , (64) 

( \ 12 

G Ts(p) = °tsl + C ts 2P* + C ts3 pz + C tsi p , (65) 

g tv(p) = °tvi + c tv2 p^ + c tv3 pi + C tVi p , (66) 

where the inclusion of the rearrangement term TP R redefines the isoscalar- 
vector coefficients: 



Cvl = 


Cvl 




C V 2 = 


f (6c„ 2 


— c s2 


C V 3 — 


\(3c v3 


- c s3 


C V 4 = 


2 (2c^4 


- C s4 



(67) 



2 small differences between p and ps at nuclear matter densities are neglected here. 
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Table 5 

Coefficients of the expansions (52-53) up to order k% in the isoscalar and isovector 
channel. The reference scale is A = 2Trf 7T ~ 0.58 GeV. 



C3 


3.01 


^3 


-2.80 


C 4 


2.67 


di 


-0.35 


C5 


2.17 


d 5 


-8.88 


C6 


13.70 


d 6 


11.50 



Table 6 

Coefficients of the expansion of the density-dependent couplings (63) - (66) in pow- 
ers of the baryon density. The input parameters in this example are those of the 
ChPT calculation [17] (63 = -3.05, a 3 = 2.16, b 5 = a 5 = 0, b 6 = -a 6 = 2.82). 



Csl — C v i 


-2.65 fm 2 


Ct si 


= Qui 


1.31 fm 2 


Cvl 


-2.65 fm 2 


C S 2 = Cv2 


1.91 fm 3 


Ctsl 


= Ctv2 


-0.25 fm 3 


C V 2 


1.36 fm 3 


Cs3 — Cv3 


1.29 fm 4 


Cts3 


= Ctv3 


-5.29 fm 4 


C V 3 


0.65 fm 4 


Cs4 = C-vAl 


2.59 fm 5 


CtsA 


= Ctvi 


6.65 fm 5 


c v4 


0.43 fm 5 



The coefficients of the expansion of the density-dependent couplings in powers 
of the baryon density (Eqs. (63) - (67)) are listed in Table 6. 
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Fig. 1. Coefficients U n in the expansion (46) of the nuclear single particle potential 
in fractional powers n/3 of the density p(r) for the case of a density p ~ p S at- 
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Fig. 2. Deviations (in percent) of the calculated binding energies (upper panel) 
and charge radii (lower panel) from the experimental values [57,58], for three rela- 
tivistic effective interactions: the non-linear meson-exchange NL3 [48] (circles), the 
density-dependent meson-exchange DD-ME1 [49] (squares), and the microscopic 
density-dependent point-coupling interaction FKVW (diamonds). The calculations 
are performed in the RMF+BCS model with empirical gaps [47] . 
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Fig. 3. Deviations (in percent) of the calculated binding energies (upper panel) 
and charge radii (lower panel) from the experimental values [57,58] for the den- 
sity-dependent point-coupling interaction FKVW with explicit inclusion of A(1232) 
excitations (diamonds) used in the present work, and for the previous version [12] 
which did not include the A(1232) degree of freedom (cirles). 
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Fig. 4. The deviations (in percent) of the calculated binding energies from the 
experimental values (upper panel) [57], and the calculated charge radii in com- 
parison with data [58], for the chain of even- A Sn isotopes. The theoretical values 
are calculated in the RHB model with the DD-ME1 [49] (triangles) and FKVW 
(dots) density-dependent effective interactions, and with the Gogny interaction in 
the pairing channel. 
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Fig. 5. The deviations (in percent) of the calculated binding energies from the 
experimental values (upper panel) [57], and the calculated charge isotope shifts in 
comparison with data [61], for the chain of even-A Pb isotopes. The charge isotope 
shifts are defined: Ar^ = r 2 ch {A) - r 2 h ( 208 Pb) and /\r 2 LD = r 2 LD (A) - r 2 LD ( 20S Pb), 
where the liquid-drop estimate is r 2 jD (A) = ^r^A 2 / 3 . The theoretical values are 
calculated in the RHB model with the DD-ME1 (triangles) [49]and FKVW (dots) 
density-dependent effective interactions, and with the Gogny interaction in the pair- 
ing channel. 
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Fig. 6. FKVW plus Gogny RHB-model predictions for the differences between the 
neutron and proton rms radii of Sn (upper panel) and Pb (lower panel) isotopes, 
in comparison with available data [62,63,64]. 
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Fig. 7. Charge form factors of 48 Ca, 90 Zr, 92 Mo, 9i Zr, 144 Sm and 208 Pb calculated 
in the relativistic point-coupling model with the FKVW density-dependent effective 
interaction (full lines), in comparison with the experimental form factors (dashed 
lines) [66]. 
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Fig. 8. The deviations (in percent) between the theoretical and experimental 
values [67] of the energy spacings between spin-orbit partner-states in doubly 
closed-shell nuclei. The calculated spin-orbit splittings correspond to the relativis- 
tic density-dependent point-coupling interaction FKVW (dots), and the relativistic 
density-dependent meson-exchange interaction DD-ME1 (diamonds) [49]. 
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Fig. 9. The deviations (in percent) of the calculated binding energies from the 
experimental values [57] of Nd, Sm, Gd, Dy, Er, Yb, Hf, W, and Os isotopes. 
The ground-state binding energies have been calculated in the RHB model with the 
FKVW parameterization, and with the Gogny interaction in the pairing channel. 
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Fig. 10. RHB model (FKVW interaction plus Gogny pairing) predictions for the 
ground-state charge radii of Nd, Sm, Gd, Dy, Er, Yb, Hf, W, and Os isotopes, in 
comparison with available data [58]. 
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Fig. 11. Comparison between the RHB model (FKVW interaction plus Gogny 
pairing) predictions for the ground-state quadrupole deformation parameters of the 
Nd, Sm, Gd, Dy, Er, Yb, Hf, W, and Os isotopes, and experimental values [69]. 
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Fig. 12. Binding energy curves of even-A Pb isotopes as functions of the quadrupole 
deformation 02- The curves correspond to RHB calculations with constrained 
quadrupole deformation. 
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